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We have examined the available latest SPS data on J/ip suppression in Pb+Pb and In+In collisions 
at 158 A GeV. Our employed model, with parameters fixed by p-p and p-nucleus collisions, gives 
excellent description of NA50 and NA60 data on centrality dependence of J /if) suppression. The 
model is then applied to predict the centrality dependence of J/ip production in Au+Au collisions 
at FAIR energy domain. A much larger suppression of J/ip is predicted. In addition the possible 
effects of a baryon rich medium on J/ip production is also investigated. 
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I. INTRODUCTION 

J/ip suppression has long been recognized as an impor- 
tant signature for the occurrence of color deconfinement 
in nuclear collisions. If QGP is produced in the collision 
zone, the cc binding potential gets shielded by Debye 
screening of colored partons leading to the reduction in 
the J/ip yield (HQ. However subsequent experimental in- 
vestigations have revealed a considerable suppression of 
the charmonium production already present in proton- 
nucleus (p+A) collisions, where QGP or more generally 
formation of any secondary medium is not expected. In 
these reactions, the produced cc pair may interact with 
the cold nuclear medium of the target nucleus, hindering 
the formation of a bound state. To quantify the nuclear 
effects, data for different target nucleus are convention- 
ally analyzed in the framework of the Glauber model |3| , 
and the suppression is expressed in terms of an effec- 
tive "absorption" cross-section crj/^- When studying 
J/ip suppression in nucleus-nucleus (A+A) collisions, a 
precise knowledge of nuclear effects is an essential pre- 
requisite to disentangle the genuine hot-medium effects. 
With this approach, both NA50 and NA60 collabora- 
tions at SPS, observed significant anomalous suppression 
of J/ip yield, at 158 A GeV in Pb+Pb 4] and In+In [l2| 
collisions respectively. Though the Pb+Pb data were 
found to be well explained by a variety of models [sT-flil] , 
with or without incorporating deconfinement scenario, 
none of them satisfactorily reproduced the NA60 data. 
The theorized origin of the additional suppression thus 
remained unsolved and debated. However in both the 
measurements, the corresponding value of GjiL was ex- 
tracted from the data collected in p+A collisions at 400 
GeV [1 31 ] - With the new measurements of charmonium 
in p+A collisions at 158 GeV where o-jj^ turned 
out to be almost twice as large as that at 400 GeV, the 
NA60 experiment reported the relative charmonium yield 
in In+In collisions to be compatible within errors with 
absorption in cold nuclear matter; an anomalous suppres- 
sion of about 25 - 30 % still remains visible in the most 
central Pb+Pb collisions. 

Till date no measurement exists on J/ip produc- 
tion in heavy-ion collisions below the top SPS energy, 
primarily due to their low production cross sections. 
The Compressed Baryonic Matter (CBM) experiment at 



FAIR [15], in GSI, Germany is planning to perform for 
the first time, a detailed study of charmonium produc- 
tion in nuclear collisions, at beam energies E\> = 10 — 40 
A GeV. The J/ip mesons produced at an early stage of 
the collisions might help in characterizing the confining 
status of the highly compressed baryonic medium, pre- 
dicted to be produced in these collisions. In our ear- 
lier work (l6j . we made an estimate of J/ip production 
cross sections in proton induced collisions at FAIR. For 
our study we employed and adapted the two component 
QCD based nuclear absorption model. Originally pro- 
posed by Qiu, Vary and Zhang (l7| . the model treats 
conventional normal nuclear suppression in an unconven- 
tional manner. Model parameters were tuned by analyz- 
ing the available data for inclusive J/ip cross sections in 
high energy proton-proton and proton-nucleus collisions. 
Aim of the present paper, is to extend our studies further 
to calculate the centrality dependence of J/ip produc- 
tion in nuclear collisions at FAIR. In addition, the pos- 
sible impact of a high baryon density secondary medium 
(confined/de-confined), expected to be produced as a re- 
sult of the collision, will also be investigated. 



II. BRIEF DESCRIPTION OF THE MODEL 

In our employed model, J/ip production in high en- 
ergy hadronic collisions, is assumed to be a factorisable 
two step process, (i) formation of cc pair, which is well 
accounted by perturbative QCD and (ii) formation of 
J/ip meson from the cc pair, which is non-perturbative 
in nature. At the leading order in a s , the partonic con- 
tributions to cc production come from two subprocesses: 
quark annihilation (qq — > cc) and gluon fusion (gg — > cc). 
With the K-factor accounting for effective higher order 
contributions, the single differential J/ip production cross 
section in collisions of hadrons h% and hi, at the center 
of mass energy yfs can be expressed as, 
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= K 



dQ 2 



da 



hih 2 



dQ 2 dx} 



x F c5 ^ J/i: (q 2 ), (1) 



where Q 2 — q 2 + 4m^ with mc being the mass of the 
charm quark and xp is the Feynman scaling variable. 
F C c^j/jji(q 2 ) is the transition probability that a cc pair 
with relative momentum square q 2 evolve into a physical 
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J/ip meson, in hadronic collisions. Different parametric 
forms have been formulated for the transition probability 
following the existing models of color neutralization. Out 
of them two functional forms namely the Gaussian form 
(F( G )(g 2 )) and power law form (F^ F \q 2 )) respectively 
bearing the essential features of the Color-Singlet [18| 
and Color-Octet [l9[ models have been found to describe 
the J/tp production cross section data in p+A collisions 
reasonably well. In p+A collisions, charmonium produc- 
tion gets affected by the prevailing cold nuclear matter 
of the target nucleus. At the initial stage, nuclear modi- 
fications of the parton densities inside the target nucleus 
affect the perturbative cc pair production cross section. 
In our analysis, leading order MSTW2008 20] set was 
used for free proton pdf and EPS09 [2l| interface for the 
ratio Ri(A, x, Q 2 ), that converts the free-proton distri- 
butions for each parton i, ff(x,Q 2 ), into nuclear ones, 
fi(%, Q 2 )- In nucleus-nucleus collisions, parton densities 
are modified both inside projectile and target nuclei. De- 
pending on the collision geometry, either the halo or the 
core of the nuclei will be mainly involved, and the re- 
sulting shadowing effects will be more important in the 
core than in the periphery. Hence the shadowing fac- 
tors have to be calculated for various centrality intervals. 
Assuming shadowing is proportional to the local nuclear 
density [H, EH, the spatial dependence is defined as: 
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Ri, p (A, x, Q 2 , s, z) = 1 + Nf(Ri(A, x,Q 2 ) - i)f^lli ) 

Po 

(2) 

where normalization is fixed to ensure that 

(1/A) J dsdzR l!f> (A,x,Q 2 ,s,z) = Ri(A, x, Q 2 ). At large 
radii, r(= \J (s 2 + z 2 )) » Ra and Ri iP — > 1, while at 
the nuclear centre, the modifications are larger than the 
average Ri. 

Once produced, the nascent cc pairs interact with nu- 
clear medium and gain relative square momentum at the 
rate of e 2 per unit path length inside the nuclear mat- 
ter. As a result, some of the cc pairs can gain enough 
momentum to cross the threshold to become open charm 
mesons, leading to the reduction in J/xj) yield compared 
to the nucleon-nucleon collisions. For both parameteri- 
zations of transition probability, the corresponding val- 
ues of e 2 , extracted from the analysis of p+A collision 
data [l6l |. exhibited non-trivial beam energy dependence. 
Lower be the beam energy, larger is the value of e 2 . In 
the present work we have used the previously found e 2 
values. 



III. ANALYSIS OF SPS DATA 

Let us now move forward to test the applicability of 
the model in describing the heavy-ion data on i /^i sup- 
pression at SPS. Fig[T] shows the variation of R'aa as a 
function of N part for In+In and Pb+Pb collisions as cal- 
culated from our model in comparison with the available 



FIG. 1: Centrality dependence of J/tp production in In+In (top) 
and Pb+Pb (bottom) collisions measured at same energy (S;, = 
158 A GeV) and kinematic domain (0 < y c .m. < 1)- Data are 
represented in terms of nuclear modification factor Raa plotted as 
a function of N par t estimating the collision centrality. Error bars 
include both statistical and systematic uncertainties. Two different 
parametric forms of the transition function are used for generating 
the theoretical curves. 



latest data [25j . The In+In data points can be reasonably 
described within errors by both Gaussian (F^- G \q 2 )) as 
well as power law (F^ p \q 2 )) forms of transition probabil- 
ity. In case of Pb+Pb collisions, (q 2 ) gives lower sup- 
pression than that observed in data. However F^ p \q 2 ) 
can fairly describe the data for all centralities and hence 
does not provide any additional room for any anomalous 
suppression mechanism to set in. For F^ G,> (q 2 ), the cor- 
responding suppression is equivalent to that obtained in 
first order approximation of Glauber theory [l6j . The 
corresponding value of e 2 was obtained by analyzing the 
recent NA60 data for p+A collisions at 158 A GeV. Thus 
it can account for the In+In data but fails to generate 
enough suppression for Pb+Pb case. On the other hand 
due to threshold effect power law form generates a much 
stronger suppression for collisions involving heavy nuclei. 
As all the model parameters are constrained from the 
p+p and p+A data, in our present calculations, no free 
parameter is required to be tuned. The observed J/ip 
suppression in Pb+Pb collisions can be fully accounted 
for by the heavy quark re-scattering in the cold nuclear 
medium, without considering further suppression in the 
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hot medium created in the later expansion stages. Earlier 
the model has also been found successful to describe the 
then available NA50 data on 3/ip suppression in Pb+Pb 
collisions Q. However in those studies shadowing cor- 
rections to nuclear parton densities were ignored and Et 
fluctuations had to be explicitly incorporated, through a 
tunable parameter, for better reproduction of the data 
at large Et- 



IV. PREDICTIONS FOR FAIR ENERGIES 

Our ultimate goal is to estimate the 3/tf> yield in nu- 
clear collisions at energies relevant to those available at 
FAIR. For this purpose, we will now use our model with 
the power law form of transition probability (F( p > (q 2 )) to 
calculate the centrality dependence of R^a f° r Au+Au 
reactions at a bombarding energy 25 A GeV. Previously 
predictions of 3/tp survival probability at this energy were 
made within transport model calculations [26j . Nuclear 
effects were incorporated through conventional Glauber 
suppression scenario. For simulating the anomalous sup- 
pression two different scenarios namely 'QGP threshold 
melting' and 'hadronic co-mover absorption' were inde- 
pendently studied. For partonic scenario, a variant of 
the geometrical threshold model @ was used with dif- 
ferent melting energy densities for different charmonium 
states. For the hadronic dissociation, inelastic collisions 
with different mesons was considered. However in those 
calculations, the magnitude of the CNM effects at FAIR 
is possibly underestimated as the value of effective ab- 
sorption cross section was taken from the p+A measure- 
ments at 400 GeV. Our present estimates predict a much 
larger nuclear suppression. Note that the degree of sup- 
pression induced by cold nuclear matter strongly depends 
on the passing time, td = 2Ra/"/, of the two colliding nu- 
clei, where Ra is the nuclear radius and 7 is the Lorentz 
contraction factor. At SPS energy (E c . m . — 17.3 GeV) 
the collision time is about 1 fm/c and the magnitude 
of nuclear effects are large. At FAIR energies, the col- 
lision time (td ~ 3 fm/c) is even much longer and the 
J/ip mesons during their evolution, will mostly encounter 
the (primary) nuclear medium rather than any secondary 
medium formed eventually due to the collision. Hence 
nuclear suppression will possibly play the most promi- 
nent role to govern the overall suppression pattern and 
we refrain ourselves to consider any probable additional 
suppression due to mesonic co-movers. However at FAIR 
energy regime, formation of highly compressed baryonic 
matter at low temperature is anticipated. Monte Carlo 
simulations [24| indicate the maximum baryon density in 
a central (b = 0) Au + Au collision at FAIR energy to 
reach as high as ps — 10po- Thus the possible imprint of 
such a high density medium on 3/tp production might be 
worth investigable. Since J /ip formation time {tj/^, ~ 0.5 
fm/c) is small compared to that required for the forma- 
tion of any secondary medium, the highly compressed 
medium will most likely encounter the color neutral phys- 
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FIG. 2: Centrality dependence of J/ip suppression at 25 A GeV 
Au+Au collisions. In addition to nuclear effects, additional sup- 
pression due to high baryon density a) confined medium (top) and 
b) deconfmed medium (bottom) are also shown. 



ical mesons rather than their precursors. The additional 
suppression induced by a confined baryon dense medium 
can be schematically expressed as: 



S p j^(b,s) = exp(- 



drp B (b,s,T) < v<tj/tp_ n >) 



In the above equation, <7,/m_jv = 6.8 mb [271 IS the aver- 
age inelastic cross section of the nucleons with the already 
formed 3/ip, v ~ 0.6 Q is 3/tp velocity and ps(b,s,r) 
is the net baryon density at proper time r at the J/t/>'s 
position. Following [Hj], the spatial dependence of the 
net baryon number density is set with the transverse 
profile of the participant density, obtained in a Glauber 
model. To and 17 respectively denotes the medium for- 
mation time and the interaction time up to which J/'i/j's 
will continue interacting with the medium. Both of them 
will depend on the path length through the nucleus and 
can be obtained from Q ■ The evolution of baryon den- 
sity with proper time r can be followed from the equa- 
tion for conservation of net baryonic current. If we ne- 
glect the transverse expansion (assuming that transverse 
expansion is slow and 3 ftp suppression occurs much be- 
fore the transverse expansion sets in), we are left with, 
toPb(t ) = tpb(t). 

The suppression pattern induced by a confined com- 
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FIG. 3: Comparison of the two different scenarios of J/tp suppres- 
sion in a compressed baryonic medium. 



pressed baryonic medium is then shown in the top panel 
of Figl^l Calculations are performed for different peak 
densities varying from po to 10pa- Higher be the density 
more violent is the suppression. If the maximum den- 
sity of the produced medium is as high as lOpo; ^aa 
approaches to zero and almost no 3/tp will survive. How- 
ever if such high density is achieved in the initial phase of 
the collision, deconfinement might set in resulting a phase 
governed by partonic degrees of freedom. In a partonic 
phase the J/'i/' will interact differently with the medium. 
The interaction potential binding the c and c together 
will be subject to Debye screening induced by the free 
color charges. To mimic the suppression pattern in a 
deconfmed plasma, we follow the geometrical threshold 
model [5], without considering the detailed microscopic 
dynamics. In this model the J/?/> suppression function, 
at an impact parameter b, can be written as: 



n p (b,s)) 



(4) 



The density n p (b, s) in the step function is proportional 
to the local energy density of the matter at position (b, s). 
In the hot and dense part of the fireball where n p is larger 
than a critical/threshold value n c , all the J/ip are ab- 
sorbed in the medium and those outside this region only 
suffer normal suppression. The threshold density n c in 
this model is a parameter, generally fixed from the data. 
However it has been observed earlier that a critical den- 
sity n c ~ 3.6 — 3.7/m -2 can reasonably describe both 
the data sets from SPS @, i] and RHIC [Hj]. n c can be 
thought of to be proportional to the threshold dissocia- 
tion energy density (ei ) required for melting of J/ip. 
If we assume a constant value of critical energy density 
(e c ~ IGeV//m 3 ), independent of baryon chemical po- 
tential p b , required for deconfinement transition, then by 
analogy the threshold dissociation energy density (e^^) 



and consequently the critical participant density, n c , can 
be assumed to be constant. The right panel of Fig|2j rep- 
resents the behavior of R J Jt for three illustrative cases 
with three different critical densities. Smaller be the crit- 
ical density, lower will be the energy density required for 
J/ip melting and more will be the suppression. We put 
an end to this section by making a comparative study for 
these above two different mechanisms of anomalous sup- 
pression. For this purpose we consider two illustrative 
cases: a) confined baryonic medium with highest possi- 
ble net baryon density {ps — 10po) and b) deconfmed 
medium with approximately constant threshold energy 
(and hence participant) density. The results are shown 
in FiglHl Two different mechanisms produce distinguish- 
ably different amount of suppressions. In a confined high 
baryon density medium, dissociation is more severe com- 
pared to that in QGP phase. Thus measurement of 3/tp 
production in nuclear collisions at FAIR might also fur- 
nish valuable information about the phase structure and 
the relevant degrees of freedom in such a high baryon 
density environment never observed before. 



V. SUMMARY 

In summary, we have estimated the J/ip production 
and its possible interactions in a high baryon density 
medium anticipated in low energy nuclear collisions at 
FAIR. Our model satisfactorily describes the J/i/j sup- 
pression data in heavy- ion collisions at SPS, with model 
parameters being fixed from p+p and p+A data. At 
FAIR, exogamous production in both the partonic and 
hadronic phase is expected to be small and the primor- 
dial production will dominate the overall J ftp yield. Con- 
sequently such measurements will offer us the golden op- 
portunity to exactly trace out the possible suppression 
pattern which will not get masked by the subsequent re- 
generation. Moreover at low energies, collision time is 
much longer and the lifetime of the produced medium 
is much shorter and nuclear effects start playing a dom- 
inant role in deciding the observed charmonium yield. 
Even at SPS, the magnitude of the nuclear effects, in 
our employed framework, are substantially large to fully 
account for the observed J/ip suppression in Pb+Pb col- 
lisions. At FAIR effects of the cold nuclear matter will 
be further amplified leading to a strong reduction of the 
J/V> yield in most central collisions. The fully formed 
J/ip mesons surviving the nuclear dissociation can sub- 
sequently interact with produced high density medium 
and undergo further suppression. The degree as well as 
the mechanism of this additional suppression depends on 
the net baryon density achieved in the collision and the 
confining status of the medium. 
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